LOCALIZED POLYNOMIAL FRAMES ON THE INTERVAL 
WITH JACOBI WEIGHTS 



PENCHO PETRUSHEV AND YUAN XU 



Abstract. As is well known the kernel of the orthogonal projector onto the 
polynomials of degree n in L^(tOQ_j3, [—1, 1]) with 10^,3(4) = (1 — t)"(l+t)'^ can 
be written in terms of Jacobi polynomials. It is shown that if the coefficients 
in this kernel are smoothed out by sampling a C°° function then the resulting 
function has almost exponential (faster than any polynomial) rate of decay 
away from the main diagonal. This result is used for the construction of 
tight polynomial frames for L'^iw^j^) with elements having almost exponential 
localization. 



1. Introduction 

A basic technique in Harmonic analysis is to represent functions or distribu- 
tions as linear combinations of functions of a particularly simple nature (building 
blocks), which form bases or frames. Meyer's wavelets [6] and the (^-transform of 
Frazier and Jawerth [3] provide such building blocks on M'^. A distinctive feature 
of Meyer's wavelets and the elements of Frazier- Jawerth is their almost exponential 
space localization and simple structure on the frequency side. This makes them an 
universal tool for decomposition of spaces of functions and distributions on R'^. 

Our primary goal in this article is to develop similar building blocks for decom- 
position of weighted spaces on [—1, 1] with weight 

and in particular, LP{wa,p)^ Hardy spaces HP{'Wc,,f3), Besov spaces, and more gen- 
eral Triebel-Lizorkin spaces. The structure of the weighted spaces on [—1, 1] is 
different and more complicated than the structure of the spaces on R due to the 
fact that there are no simple translation or dilation operators in these spaces. This 
creates a great deal of complications. 

The Jacobi polynomials {Pn"'^^}^^^ provide a basic vehicle for representation 
and analysis of functions or distributions in the weighted spaces on [—1,1] with 
weight Wa^i3. We let Cq_^ denote the normalization constant of Wa.p, i.e. := 

J-i Wa,f3{t)dt. The Jacobi polynomials {Pn"''^^}'^=o are orthogonal with respect to 
Wa,f3, namely [8], 

Date: February 2, 2008. 

1991 Mathematics Subject Classification. 42A38, 42B08, 42B15. 
Key words and phrases. Localized polynomial, Jacobi weight, frames. 

The first author has been supported by NSF Grant DMS-0200665 and the second author by 
NSF Grant DMS-0201669. 



2 



PENCHO PETRUSHEV AND YUAN XU 



where 

r(a + /3 + 2) r{n + a + l)T{n + (3+l) 



r(a + l)r(/3 +l){2n + a + P + l)r(n + l)r(n + a + /3 + 1) 
For / e L^(?«Q_^) the Fourier expansion in Jacobi polynomials is 

n=0 -^-1 

The nth partial sum of the expansion can be written in terms of the reproducing 
kernel K^'^\x,y) as 

n „1 

Sn{f;x) =Y,dAmf'^^r^P^'''^\x) = c„,^ / f{t)Kl^'P\x,t)w^,fi{t)dt, 
where the kernel is given by 

(1.1) K^-''\x,y) = {^-''Y pt-''i-)Pt'''(y)- 

One of our main results in this article asserts that any polynomial in two variables 

of the form 

°° / ■ \ —1 

(1.2) K'^ix,y) = U) {h^^^f j>(".«(,)^(-/^)(,), 

3=0 ^ ^ 

where a e C°°[0, oo) with supp a C [1/2, 2] (or more generally suppa C [c, d], c > 0) 
has almost exponential localization around the main diagonal y = a; of [— 1, 1]^. To 

be more precise, let 

Wa,/3(n; x) := (1 - X + n-2)"+i/2(i + 3. + „-2)/3+i/2^ -1 < a; < 1, 
then for a, /3 > —1/2 and any A; > 1, there is a constant Cfc > such that 

(1.3) |Lg'^(cos6l,cos^)| < Cfc^ = 

\Jwa,ji{n\co&Q)^Jwa,ii{n\ cos(^)(l + n\Q - 

for < 0, (/) < TT. This result is a far reaching extension of the recent discovery in [7] 
that L^''^{x, 1) (with A a half integer) has almost exponential localization around 
X = 1, which is utilized in [7] for the construction of frames on the n dimensional 
sphere. The above result shows that if one smooths out the coefficients of the 
kernel in (1.1) by sampling a C°° function then the resulting function decays away 
from the main diagonal at almost exponential (faster than any polynomial) rate. 
This fact, which is well known for the trigonometric system, has apparently been 
overlooked. We believe that it will play an important role in various problems for 
weighted spaces. 

The polynomials L'^'^ give us a handy tool for constructing tight frames in 
L^{wa,i3) and other weighted spaces (with weight Wa.jj) on [—1, 1]. Our construc- 
tion of frames utilizes a semi-discrete Calderon type decomposition coupled with 
discretization using the Gaussian quadrature formula (see §3). A similar scheme is 
used in [7] for the construction of frames on the sphere. 

Let us denote hy ip^, ^ X , the constructed frame elements (§3), where X — 
Ujlo'^j is a multilevel index set consisting of the localization points (poles) of the 
V'^'s. Then our frame is defined by 
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We show that every function / e ^^(wq,/?) has the representation 

i.e. is a tight frame in L^{wa,0)- 

A distinctive property of the jth level frame elements tp^ e Xj) is their local- 
ization: 

(1.4) \Mcose)\ < cfe nsnloiia ^ 

^ywa,0{2^ ; cos + 2^9 — arccos^D* 

which is analogous to tlic locaHzation of Meyer's wavelets. It is worthwhile to 
point out that the factor y'wa,/? (2-' ; cos 0) in the denominator in (1.4) as well as 
the corresponding terms in (1.3) reflect the expected "effect" at the end points of 
[—1,1] and play a critical role in our development. Notice that this term is present 
even in the case when Wa,/3{t) = 1 (a = /? = 0). 

The superb localization of the frame elements {V'j} prompted us to term them 
needlets. This along with the semi-orthogonal structure of 5* and the increasing 
number of vanishing moments of the tp^ 's enables one to utilize the needlet system ^ 
for decomposition of spaces other than L'^{wa,i3) such as Besov and Triebel-Lizorkin 
spaces. Wc; shall report on results of this kind in a follow-up paper. 

We do not discuss here the literature on polynomial bases and frames on the 
interval because to the best of our knowledge the elements of the existing bases 
or frames do not have the localization of the needlets. It is an open problem to 
construct bases in L'^{wa,fi) with basis elements having localization similar to (1.4). 

This article is organized as follows. In §2, we establish the localization (1.3) of 
the polynomials L'^^^ from (1.2), where in §2.1 we consider the the particular case 
(/) = (j/ = 1) , while in §2.2 we prove (1.3) in general. In §3, we construct our 
polynomial frames and establish their main properties. 

Throughout the article positive constants are denoted by c, Ci , . . . ; unless spec- 
ified, their values may vary at every occurrence. The notation A B means 
ciA < B < C2A, and A := B ot B =: A stands for "A is by definition equal to B" . 

2. Localized polynomials in terms of Jacobi polynomials 

In this section we establish the localization properties of L!^'^{x,y) from (1.2) 
depending on the smoothness of a. We first prove (1.3) when ^ = {y = 1) and 
then consider the general case. 

2.1. The localization of L"'^{x, y) in the case y = 1. 

Theorem 2.1. Let a G C'^[0, 00) with k > 1 and suppa C [1/2,2]. Assume that 
a> (3 > —1/2 and define 

(2.1) L^''^(cos^) = E« U ('^^''O" ^^"'''(l)^^"'''(cos^). 

Then there exists a constant > depending only on k, a, (3, and a such that 

(2.2) K-"(cos^)|<c. ^^^^^^,^^_^ , 0<^<7r. 

The dependence of Cfe on a is of the form Cfe = c(a, /3, k) maxi<y<fc jja^"^) jj . 
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Proof. Since P^"'''^ (1) = (a + l)„/n! = T{n + a + l)/[r(a + l)r(n + 1)] [8, (4.1.1)], 
it is easy to verify that 

(2.3) ^n'(^<^^^)-r^^Tl3h 

X ±a (I) (2. + a + /3 + l)r(,+c. + /? + l)^(.,,) 

^ U; r(j + /3 + i) ^ ^ ' 

Notice that since suppa C [1/2,2], the sum in (2.1) (and (2.3)) is finite; in fact it 
is over n/2 < j < 2n. 

We first prove (2.2) for < ^ < tt/u. Using the fact that [8, Theorem 7.32.1] 
ll^'i"''^^l|L~[-i,i] ~ cj" and r(j + a)/T{j + 1) ~ j"'^ as j oo, it follows that 

2n 



]L«''5(C0S^)1 < C^i^a+l < ^^2a+2^ 
3=0 

which yields (2.2). 

Let w/n < 6 < TT. We will use the identity [4]: 



(a0) =(ia,/3(l-cos^)-" / COS [m + (a + /3 + l)/2](/) 

-Pm (1) •'O 

{cos (j) — COS 6)"^^/'^ fa + f3a — (3 1 cos0 — cos(/)\ 
(l + cos0)("+«/2 1,^' " + 2 ' 1-cos^^ j 

for < 9 < TT and a > —1/2, where 2F1 denotes the hypergeomctric function 
and da,0 = 2(°+^+i)/2r(a + l)/(V7Fr(a + 1/2)). Recall the identity [8, (4.1.3)] 
Pm'^\x) — {—l)"^Pm'°'\—x). Now changing first the variables i-^- tt — in the 
integral, and then setting 9 n — 9 and interchanging the places of a and (3, we 
obtain 

(2.4) "",3^ ^ = rf/3,c«(l + cos^)-'^ / cos [m<^ - (a + /3 + l)(7r - <^)/2] 

-Pm ^1) 

(cos^ — cos^)^"-^/^ /a + /3/3 — a 1 cos(^ — cos6'\ 

(l-COS<^)(«+/3)/2 ^ + 2' 1-COS^ j 

Since pf'''\l) = T{j + p+ 1)/[T{(3 + l)r{j + 1)], we have 

„(/3,a). . r(/3 + i) r(j + a + /3 + i) ^ 1 r(j + a + /3 + i) 

^' ^ ^r(a + /3 + 2) r(j+/3 + l) r(a + /3 + 2) r(j + 1) 
Using this and combining (2.3) with (2.4), we get 

L^'^icos9) = dl^f,{l + cos9)-P r [cos(A0 - ATr)^^ (<^) - sin(A</. - A7r)^f (0)] 

J9 

(2.5) 

(cos0 - cos(/>)^~^/^ fa + P (3 - a 1 cos0-cos^\ 
^ (l-cos<^)(«+/3)/2 (,^'^'^+2' 1-cos^ 

where := 2d/3,a/r(a + /3 + 1) and 

,cos..^ ._ V ^ A (J + A)r(i + 2A) 
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with A := (a + + l)/2. 

The idea of the' proof is to derive (2.2) from tlic; analogous loealization of the 
trigonometric polynomials and A^'^. This in turn will follow by employing the 
fact that any trigonometric polynomial An{9) = J2jLi aje^^^ has an excellent local- 
ization around zero whenever the coefficients aj come from sampling of a smooth 
compactly supported function. 

Define 

^ ^^ft + A)rft + 2A) ^^^^^ 

^' r(t + i) 

Notice that 2A = a + /3 + l>0, since a> j3> -1/2. 

With the next lemma we establish the smoothness and decay as f — > oo of G and 
its derivatives. 

Lemma 2.2. The function G is analytic on [0, oo) and for any k>Q 
(2.6) \G^''\t)\ < Ckt^^-'' , t>l, 

where Ck depends only on k and A. 

Proof. We first show that for 6 > a > 0, cr := a + 1 — 6, and t > 0, 

(2-7) iS^ = (-im(a)^^fe^^(l + c:^r^)) ast-.oo, 
^ ' dt''T{t + b) ^ ' °^ ' T{b - a) t^-«-+^ \ ^ \ jj 

where B'^{u) (here n = 0) are the generalized Bernoulli polynomials defined in [1, 
(C.4.2)]. We will slightly modify the argument from [1]. 
The relation between Gamma and Beta functions gives 



1 



r(6 - a) Jo 



It is readily seen that the last integral is an analytic function for f > and it can 
be differentiated with respect to t as follows 



d>' T{t + a) _ , ^.st^k^-san „-M6-a-l, 

dt'' T{t + b) T{b - 



a) Jo 



Now exactly as in [1, C.4, p. 615] 



(1 - e-^y n! 

and by Watson's lemma, we obtain the asymptotic expansion 

Tjt + a) ^ {-l)k+n r{b-a + k + n) 1 
dtT{t + b)^^^ n! ^ r{b-a) 

which implies (2.7). 



oo, 
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Let a :— 2X — [2AJ , where [x\ denotes the largest integer less than x. Then 
< a < 1. Using that r(l + x) = xr{x), we have 



G{t) = 



(t+l)(t + A)r(t + 2A) 

r(i + 2) 



{t + m + X)it + 2X-l)...{t + a)S^ = p{t)^^' + 



'T{t + 2) '"■"'T{t + 2)'^ 
where p{t) is a polynomial of degree [2AJ + 2. Now by (2.7), 

T{t + a) 



and hence 



dt" T{t + 2) 



< 



i2- 



— — , t>l, iy = 0,1,. 



v=0 

k 



dt^-" T{t + 2) 



L/=0 



y j -^2 — a-\-k—u 



< ct 



2X-k 



Thus (2.6) is established. 



□ 



Now relying on the smoothness and localization properties of a and G, we derive 
the desired localization of A'^^ and A^'^. 

Lemma 2.3. If A^ = A^°^ or A„ = then 

(2-8) K(^)|<c,^3^^, 1^1 <7r, 

where Ck is of the form Ck = c{a, (3, k) maxi<y<fe \\a'^'^'' H^i with c(a, j3, k) depending 
only on a, /3, and k. 

Proof. Let n > 2 and set 

A„(e) :=^a(|i|/n)s(i)G(U|)e'^^ 

where s{t) = 1 or s{t) = signi. Evidently, A„(^) = 2A^^%e) if s(i) = 1 and 
A„(6') = 2iA^^''{9) if s(<) = signt. 

Set $„(t) := a(|t|/n)s(t)G(|i|) and define i>„(i) := ^ 4$„(0e^*«dC- Clearly, 
$„(i) e C''(M) and supp$„ C [-2n, -n/2] U [n/2,2n]. Consequently, 
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and using (2.6), 



j=0 



2^ Q ^ ^ |a«(«)G('=-^)(nw)|du 



<cn2^-'^+i max ||a« l^i = Cfen^^-'^+i. 

0<j<fe 

Therefore, |$n(i)| < Cfcn^^~'^+^|f|~'°. As above (with fc = 0) but easier, we obtain 
|^ra(i)| < cn^^+^. Combining these, we get 

(2.9) |$„(i)|<c,^^--^, tern. 

Applying the Poisson summation formula 

27r^/(27ri) = ^7(j) 

iez jez 

to the function f{t) := $„(t + 6*) (then f{t) = $„(t)e"''), wc get 

A„(0) = ^$„(j)e^J"^ = 27r^$„(27rj +0). 
jez jez 

Now using (2.9), we obtain for |^| < n, 



jez jez 



(l + n|e|)'= ^ (1 + n(2j - l)7r)'= 



< c- 



„2A+1 



'(l + nl^l)'^ 

and estimate (2.8) follows. □ 

We are now in a position to complete the proof of estimate (2.2) in the case 
n/n < 6 < TT. Since < (cos^ — cos(f))/{l — cos6) < 1 for ^ < < tt, the absolute 
value of the 2-F1 term under the integral in (2.5) is bounded by a constant. Hence, 
using (2.5) and Lemma 2.3, we infer 

(2.0) iL-(c„«)i < c(i + c„.,-^ j- ;;";;'°r-w;lvr * 
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In the following we will use the identities cos 0— cos = 2 sin sin 1 — cos (}> = 
2 sin^ |, and 1 + cos </> = 2 cos^ |. 

Consider first the case ■K/n<6< 7r/2. Then for ^ < ^ < tt, 

2\/2,^ . e + (i) + (j) , (p-0 . (j)-e (p-e 

-^{e + 6)<sin r < _Lr and <sin^^ < -^^ . 

^■K ^ ^' - 2-2 -K - 2-2 

Using these in (2.10), we get 



|L^'^(cos^)| < cn«+^+2 /" 
Je 



dcj) 



(02 _ 02-^0-1/2 



„2a+2 .oo .2 _ l^0-l/2 ^2a+2 

< C^^rrrr ^ / vrrh^ du < c- 

Ji 



In the case 7r/2 < 6 < n, estimate (2.10) evidently gives 

i.cS/ nM n«+^+2 . r (cose-cos<^)'5-i/2 

L"''5(cos6') < c- ^(1 + cos0)-^ / ^^^-^ — 

' " ^ - (1 + 71^)*=^ ^ Jg (l-cos0)«+/5 ^ 

Let 6*' = TT - 61. Then < 6*' < 7r/2. Hence, for < </> < i9', (i9' ± 0)/7r < 
sin[(0' ± <j))/2] < {9' ± (f))/2. Consequently, substituting <f> i-^ n — <j) gives 

r'^ (cos 9 - cos ^!))/3-i/2 /-e' (cos - cos 9'f-^/^ , , 

-'-^ — ' -dcp 



(1 - cos (/))«+/' Jq (1 + C0S(^)"+'' 

Jo 

Since 1 + cos6i = 1 - cos6i' = 2sin^(6i72) ~ {9'f, this shows that 

^a+0+2 n^«+2 

1^" ^ ^(TTW ^ ^ (l + n^)^+^-/^ 
as 7r/2 < 6 < TT. The proof of Theorem 2.1 is complete. □ 

2.2. The localization of L"'^{x,y) in the general case. Recall the definition 
of L^''^(a;,y) from (1.2): 

LZ'^ix,y) ■.= ±aU) (/.;-'^))"^P,^-^)(.)P,(-''^(y). 

In this subsection we estimate the localization of L!^'^{x, y) around the main diag- 
onal ?/ = X of [— 1, 1]^, which depends on the smoothness of a. To this end we will 
need the quantity 

(2.11) «;„,^(n;a;) := (l-a; + n-2)«+V2(i + a; + „-2)/3+i/2^ -1 < a; < 1. 
Notice that since sin| sin^ on [0,27r/3] and cos | ^ sin^ on [7r/3,7r], then 

(2.12) Wo,,0{n;cos9) ~ (sin6l + n-^)^"+\ < 61 < 27r/3, 
and 

(2.13) Wo,,0{n;cos9)r^{sm.9 + n-'^f^+'^, 7r/3 < 6» < tt. 
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Theorem 2.4. Let a, (3 > -1/2 and let a e C'=[0, oo) with fc > 2a + 2/3 + 3 and 
suppa C [1/2,2]. Then there is a constant Ck > depending only on k, a, /?, and 

a such that for < 9, 4> < tt 

(2.14) \LZ^''{cos9,cos^)\ < cu T ^ m ^h-^ ' 

where u = k — 2a — 2p — 3. Here the dependence of Ck on a is of the form Ck = 

c(a, /3, k) maxi<^<fc Ha^"^) || . 

Proof. We need the product formula of Jacobi polynomials [5]: For a>/3>— 1/2, 

Pi"'^'(x)Pi"'^'(y) 



P("''^)(l) ~'^Jo Jo 

where 

i(x, y, r,il)) = \{l + x){l + y) + \{l - x){l - y^"^ + ry/l - x^^l - cos^" - 1, 
the integral is against 

dm(r, V) = (1 - r2)«-'5-V2'3+i(sin^)2/'dr#, 
and the constant Ca,/3 is selected so that 

Ca,/3 / / ldm{r,i!) = 1. 



= Ca,0 



Pi^'^\t{x,y,r,i,))dm{r,i^), 



/O JO 

Once (2.14) is established for a > /? > — 1/2, the case (3 > a> —1/2 will follow from 
the relation P^^'^^^cos 6») = (-l)"P^'''"^(cos(7r - 9)). In the case a = /3 = A - 1/2 
we use the product formula of Gegenbauer polynomials [2, Vol I, Sec. 3.15.1, (20)] 

a^(cos6')C<^(cos?!>) ^ , , 

-J-^ — ' -1^ — = cx I Cf{cos9cos(l) + usm9sm4,){l-u^)^-^du, 

Cj (1) J_i 

where c\ is the normalization constant of the weight function {1 — u^)^~^ . Since 
the case a = /? is much easier, we shall only give the proof when a > p. 

The product formula allows us to write L^'^{cos9,cos(p) in terms of L^'^{t) 
defined in Theorem 2.1, namely, 

(2.15) Ll'^{x,y) = c^,J f LZ'^{t{x,y,r,^))dm{r,i,). 

Jo Jo 

For t = cos 9 with < ^ < tt, we have 9 ^ sin^/2 ^/l — t, and consequently 
the estimate of \L"'^{t)\ from Theorem 2.1 can be rewritten as 

We will apply this estimate with t = t{x,y,r,xp). Let x = cos 9 and y = cos(/>, 
< 61, < TT. Then 

1 — t{x, y, r, ip) =1 — cos 9 cos cp — sm9 sin (f) 

(2.17) + 5(1 - cos6')(l - cos(/))(l - r'^)+sm9sm(f){l -rcosip) 

n 9 — d) 2^20/ 2\ /\ 

=2sin — ^ |-2sin -^sin — (1 — r ) + siny sin(j!)(l — rcos'^), 

yielding 

1 - t{x,y,r,i^) > 2sin2 ^ ^{9- c^f. 
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Then by (2.15)-(2.16), wc get 



(2.18) |i«./5(cos^,cos<^)|<^^^^^^-^^^^, 0<^,<^<7r. 

Our next step is to show that this estimate can be substantially improved away 
from the end points of the interval. 

Case (a): 0<9< 27r/3 and < < 7r/2. By (2.17), we have 

1 — t{x,y,r,tp) > c{\9 — 0p + sin sin 0(1 — rcostpj). 
Then by (2.15)-(2.16), 

Jo Jo (1 + n[|^ — + sm6'sin0(l — r cos'^)J^/^j 

wherea = k + a-l3. Set F{u) := 1/(1 + n(|6i- + sin6isin 0(1 - u))V2)<^. in the 
following we evaluate the double integral by changing the variables several times: 

r [ F(rcos^)dm(r,V)= / / F{rs){l - r'^)"-^-^r'^'^+'^{l - sy-^^^dsdr 
Jo Jo Jo J-1 

= j'^ F{u)r{l - r2)«-/3-i(,.2 _ uy-^/^dudr 
= t F(u) f r(l - r^r-('~\r^ - u'f-^l^drdu 

J-l J\u\ 



; F{u){l - u'^^-^'^dv 



where c = (1/2) vP-^''^{l - v^'^'^dv (here v = {r^ - v?)/{l - v?) is the last 
substitution). This gives 

We now split up the last integral into two integrals: one over [—1,0] and the other 
over [0, 1]. For the first integral, 1 — u > 1 and hence 



fO 

J F(u)(l-u2)«-V2^y < 



< 



(1 + n[\9 - 0|2 + sin 61 sin 0] 1/2)0 
c 



(n2 sin6'sin(/))"+i/2(l + n\9 - 0|)fe-"-/'-i ' 
To estimate the integral over [0, 1] we apply the substitution t = 1 — u and obtain 

(1 - w)"-i/2dy 



/ F{u){l-u^Y-^'^du<c f 
Jo Jo 

I 



(1 + n[\9 - (/-|2 + sin6isin0(l - u)]^/^y 

sin sin ^a— 1/2^^ 



- (n2 sin6isin0)«+V2 (i + [n'^{0 _ 0)2 + t]i/2)a 
^ c .00 ^a-1/2^^ 

- (n2sin6'sin0)"+i/2(i + n|6'-0|)<^-2«-2 Jo (l + t)«+i 
c 

(n2 sin6isin0)«+V2(i + n\e - 0|)fe-«-/3-2 ' 



< 
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Putting these estimates together, we get 
(2.19) |L^''^(cos6l,cos<^)| < 



(sin e sin 0)«+i/2(l + n|6' - (f)])''-"-^-'^ ' 



Our last step here is to show that estimates (2.18)-(2.19) imply (2.14). Observe 
first that since < 0,(f> < 27r/3, we have by (2.12), 

(2.20) Wa,/3{n] cos e)wa,/3in; cos (I)) ~ {sine + n-^f"+\sm4, + n-^f"+^ 

~(0 + n-i)2"+i(</, + n-i)2"+i. 

We need the following simple inequality 

(2.21) (6i + n-^)((/) + n-i) < 3(6'^ + n-2)(l +n|6»- (/)|), 6»,0>O. 

To prove this we assume that (j) > 9 and define 7 > 1 from (f> = j0. Suppose 
7 > 3. Then (2.21) wiU follow if we show that {0 + (j))n-^ < 3n-'^n\0 - <p\, which is 
equivalent to7 + l<3(7— 1). But this holds since 7 > 3. 

Let now 1 < 7 < 3. Then it suffices to show that {0 + (t>)n-^ < 2{0(j) + n'^). In 
turn this inequality holds if Wn~^ < 2{0^ + n~^), which is apparently true. Thus 

(2.21) is established. 

It is now readily seen that (2.18)-(2.19) yield (2.14) by using (2.20)-(2.21). 

Case (b): tt/S < 9 < tt and 7r/2 < (f) < tt. In this case, we have by (2.17) 

(2.22) 1 - t{x, y, r, V) > c(|6» - (/>|^ + (1 - r) + r(l - cos V') sin 9 sin (j)) 
upon writing 1 — r cos tjj = 1 — r + r(l — cos'^). In particular, 

l-t{x,y,r,i;)>c{\e-<j>f + {l-r)), 
which shows, using (2.15)-(2.16) , that 

|L«'''(COS0,COS<^)| < Cn2"+2 /' (1 _ r2-)a-/3-1^2/3+l^ 

' " ^ ' Jo (l + n[|e_0|2 + (l_r)]l/2)-' 

with a := k + a — /3 as before. Applying the substitution n^(l — r) = s, we get 

(2.23) IL^^f^ {cos 9, cos cl))\ < cn^"+^ t (1 - rY-^-^dr 

^ ^ ' " ^ io (l + n(|e-<^|2 + (l-r))V2)'^ 



< c- 



i (1 



/o {l + {n^\9-4>\+sY/^Y 

^2/3+2 .00 s<^-li-lds 



(1 + n\0 - </,|)<--2(a-/3)-l/2 (1 _^ s)a-/3+l/4 

< C- 



„2/3+2 



'(l + n|6i-?i|)'=-«+/'-V2- 
Furthermore, using (2.15)-(2.16) and (2.22), we obtain 
|i°''^(cos6i,cos?!>)| 



(1 - r-f-'^-^-r-^+^[^v)-\hrdr 



Jo (l + n[|6»-(/)|2 + (l-r)+r(l-cosV')sin6»sin<?!)]i/2)'^ 
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The inner integral can be estimated by changing the variables. Set A := \9 — 0p 
1 - r. Then 



Jo "io (1 



(1 - sf-^-(ls 



+ n[A + r(l - s) sin (9 sin 0] 1/2)'^ 

-j^ ^n^rsinesin<^ ^,9-1/2^^ 



< 



< 



(nV sin6'sin0)/3+i/2 (i + (ji^A + t)i/2)<T 

1 r t^-^'^dt 



(n2rsin6lsin^!>)/3+i/2(i+nAi/2)<^-2/3-2 + 

c 



Jo 



[n'^r sin B sin (/))/5+i/2 (l + „^i/2)fe+a-3/3-2 ^ 
and since 1 — s > 1 for — 1 < s < 0, 

(1 - sf-^l^ds 



< 



7_i (1 + n{A + r(l - .s)sin6lsinc6)°' " (1 + n(A + r sin6'sin(/))V2)'T 
c 

- (nV sin e sin 0)/3+l/2 (l + „^l/2)fe+a-3/3-2 • 

Combining these, we get 

„2a+2 |-1 _ ^2^a-/3-l^/3+l/2^^ 



|L^'''(C0S6',C0S(/))| < 



(n2sin6lsin0)/3+V2 (i + _ 0|2 + (i _ ^))i/2)^ ' 



where t := /c + a — 3/3— 2. The last integral (with a in place of r) has already been 
estimated in (2.23). Consequently, 



|L^'^(cos6l,cos(;i)| < 



C„2/3+2 



(n2 sin 6* sin (/))/'+i/2(i + n|6l - 0|)fe-"-/3-5/2 ' 



As in Case (a), this estimate along with (2.23) implies (2.14), using (2.13) and 
(2.21). 

Case (c): Either < 6* < 7r/3 and 7r/2 < ^ < tt or 27r/3 <^<7rand0<(^< tt/2. 
In this case, we have \6 — (l)\> 7r/6, so that 

1 - t{x, y, r,'il})>c\e-(l)f>c> 0. 

Hence, by (2.15)-(2.16), 

I L"''^ (cos 6*, cos 0)1 < ,^ ^"x,, 3- < 7Z ^T" 3— T) 

which yields (2.14) since ti;a,/3(n; cos ^)wa,/3(n; cos 0) < c. The proof of Theorem 2.4 
is complete. □ 

The estimate of \L"'^{x, y)\ from Theorem 2.4 allows us to control its integral. 
Proposition 2.5. Let k > 3a + 3(3 + 5. Then for 1 < p < oo, 



(2.24) j'^\LZ'^{x,y)\Pw„Ay)dy< 



p—i ' 



-1< a; < 1. 
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Proof. We will prove (2.24) only in the case < x < 1, since the case — 1 < a; < 
is the same. Let x =: cosO and denote I[a,b] '■— \Ln'^{x,y)\P'Wa,0{y)dy. We first 
estimate /[_i/2,i]- Using Theorem 2.4 and (2.12) and substituting y = cos^, we get 

r2-/3 (sin0)2"+id0 



np r 



' {wa,p{n;x))P/^ Jq (sin(/> + n-i)("+i/2)p(i + - (f>\)''P 

with a = k — 2a — 2(3 — 3. Denoting the last integral by J, we use the relation 
sin ^ ~ (/) for < (/) < 27r/3 and apply the substitution u = n6 to obtain 



J <c 



n2«+i Jo (1 + n0)(a+i/2)(p-2) (1 + n\e - (t)\yp 



^C-„2a+2 (l+'u)("+l/2)(p-2)(l + |u_„5l|)ap- 

If p > 2, then we split the last integral into two integrals: one over [0, n9/2] and 
the other over [n6/2, 2n7r/3]. We have 
rne/2 ^ .ne/2 



Jo ~ (l + n6l)<^P Jo (1 + m)(«+V2)(p-2) - (l + n6l)'^P-i 



and 



2n,r/3 ^ .2n^/3 

< ,., / : TT^ < 



g/2 (1 + n0)("+i/2)(p-2) (1 + |m - n6l|)'^P " (1 + n6')(«+ V2)(p-2) ' 

using that the last integral is bounded by J^^{^ + \u\)'''^Pdu < c. 

The assumption k > 3a + 3(3 + 5 yields ap — 1 > a{p — 2) and hence 

^(a+l/2)p ^ 

J <c , < 



„2a+2(l ^„5/)(a+l/2)(p-2) - ^ „-l)(a+l/2)(p-2) 

c 

< 



n{wa,j3{n; cos 0))p/2-i 
If p < 2, then we have 

^(a+l/2)p /.2nV3 + y)(a+l/2)(2-p)^^ 



J <C 



„2a+2 (1 + |u _ „5»|)<Tp 

„(a+l/2)p .n(2V3-9) (1 + + n^)(a+l/2)(2-p)^„ 

J— 1 



Using the inequality {A + By < 2"{A" + B") we conclude that 

^(a+l/2)p 



J<C 



n 



2a+2 



^ (^^^(a+l/2)(2-p) /""^ 



, (1 + |i;|)'^P-("+l/2)(2-p) ^ ^ (1 + |^|)ap 

< c "' , ^ (1 + n^)("+l/2)(2-p) < 



^(a+l/2)p 



'^(M'a,/3(n;COS0))P/^ ^ 



which is the same as in the case of p > 2. Here we have used that k > 3a + 3(3 + 5, 

which yields ap - {a + 1/2) (2 - p) > 1. 

Putting these estimates together, we conclude that 

nP J nP~^ 
(2.25) 7[_v2,il < ^(^^^^(^.^))p/2 ^^(^„,,(„;^))P-i- 
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To estimate /[- 1,-1/2] we again apply the substitution y — cos0. Notice that 
since < 6* < 7r/2 and 27r/3 < (/> < tt, we have \9 — (j)\ > 7r/6. Employing 
Theorem 2.4 and (2.13), we obtain as above 

nP n^P+U^)p r (sin(/.)^/3+i 

[-''-'/'I - ''K,/3(n;x))f/2 ^ap y^^^g (l + nsin<^)(/3+V2)p'^'^ 

< -, ZT—u: ,„ , ., 7T < C- 



'(wa,/3(n;a;))p/2 n<^P „(/3+i/2)(p-i) - {w^^^{n;x))P-^' 

using the trivial inequalities n~'^'^°'~^^^ < Wa^i3{n]COs9) < c for < < n/2 and 
the fact that k> 3a + 3/3 + 5. The last estimate coupled with (2.25) yields (2.24). 
The proof of Proposition 2.5 is complete. □ 

3. Tight polynomial frames in L'^{wa,f3) 

We now turn to the construction of polynomial frames. In addition to the lo- 
calized polynomials from §2, we shall need the Gaussian quadrature formula. Let 

= cos9j, 1 < j <n, denote the zeros of the Jacobi polynomial Pn"'^\t), ordered 
so that 

=: 6*0 < 6*1 < • • • < 6'n < 9n+i := n 
Let A„(t) be the Christoffel function and 6^ = A„(^j,). It is known that 

^i^+i — 9v and hence 61, i'n~^ (1 < < n); 

and also 

Here the constants of equivalence depend only on a, (3 (cf. [9, p. 282]). These 
quantities appear in the well-known Gaussian quadrature: 

Proposition 3.1. For each n> 1, the quadrature 

/I n 
f{t)wc.,p{t)dt^y^Kf{^^) 
-1 .=1 

is exact for all polynomials of degree 2n — 1. 

We shall utilize Proposition 3.1 with n = 2^ . Let us denote by ^j,^ and bj^i, 
[v = 1, 2, . . . , 2-') the knots and the coefficients of the Gaussian quadrature when 
n = 2^ , and set 

Xj := : u = 1,2, . . .,2^}, and b^:=bj^i, if^^^j,^. 

Then the Gaussian quadrature 

(3.2) c„,/3 /" /(t)«;„,/3(t)c;t~ ^ 6j/(0 

is exact for all polynomials of degree at most 2^+^ — 1. 
Let a satisfy the conditions: 

(3.3) aeC°°{R), o>0, suppa C [1/2, 2], 

(3.4) a{t) > 0, if t e [3/5, 5/3], 

(3.5) a^{t) + a^{2t) = l, if ie [1/2,1]. 
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It follows from these conditions that 

oc 

(3.6) ^a2(2-''t) = 1, te[l,oo). 

It is easy to construct functions that satisfy properties (3.3)-(3.5). Indeed, a 

standard construction (e.g. in wavelets) shows that there exists a function g with 
the properties: g e C°°(R), suppg = [-1,1], git) > on (-1,1), g{-t) = g{t), 
5(0) = 1 and \g{t)\'^ + \g{t + = 1 on [-1,0]. Then the function a(t) := g{log2t) 
has the required properties. 

We introduce the orthonormal Jacobi polynomials P^'^\ which are given by 

PI^''^Hx) := (/l(f'^^))-l/2p^(a,/3)(3.)^ 

The sequence Q is an orthonormal basis in L^{wa,f3j [—1)1]) with inner 

product 

{f,g)~Ca,f3 j .fit)g{t)Wa,l3it)dt. 

Assuming that a satisfies conditions (3.3)-(3.5), we write Lq'^ := 1 and 

oo 

(3.7) Lf^ix, y):=J2a(^) P(«'« (a;)p(«'« (y), j = 1, 2, . . . . 

We now define our frame elements (needlets) by 

(3.8) ,p^ix):=^^-Lff'ix,0 for ^ e A",, j = 0, 1, . . . , 

where Xj is the set of knots and the b^'s are the weights in the Gaussian quadra- 
ture (3.2). 

Denote X := Uj^^Xj, where every two points ^,r] £ X from different levels 
Xj / Xk are considered to be different elements of X even if they coincide. We use 
X as an index set in the definition of the needlet system 

which as will be shown in Theorem 3.3 below is a tight frame in L^(u>q,^). 

We next show that each tp^ has nearly exponential (faster than any polynomial) 
rate of decay away from ^. This property of needlets is critical for their implemen- 
tation to decomposition of spaces other than L'^{wa,0)- 

Theorem 3.2. For any k > 1 there exists a constant Ck depending only on k, a, (3, 
and a such that for £, G Xj, j = 1,2, . . . , 

2i/2 

(3.9) \M^os9)\ < ck^===——- — , 6 e [0,7r], 

yjwa,i3{2^ ',cosO)[l + 2-J|6' — arccos^l)'' 

and, consequently, for v — 1,2, ... ,2^ , 

Proof. Estimates (3.9)-(3.10) follow immediately from the estimate in Theorem 2.4, 

taking into account the simple fact that 1 + 2-' |6' ± c2~^ | - 1 + 2-''6' if < 6* < tt 
and the fact that ~ 2~^Wa,i3{2~^;i). Since by assumption a is in C°°, estimates 
(3.9)-(3.10) holds for every A;'> 1. □ 
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Remark. It is worthwhile to mention that the estimate (3.9) or (3.10) (also Propo- 
sition 2.5 with p = 2) immediately yields lU^Cwc p) < c. 

Our next theorem shows that is a tight frame in L'^{wa,0)- 

Theorem 3.3. If f e L'^{wa,p), then 

oo 2' 

(3.11) / = EE(/'^«.-.^)^«.-.^ = E</'^«)^« ^'(^"'Z') 

and 

(3.12) ll/IU^(«,„,.)=(El</'^«)l')'^'- 

Proof. Suppose / e L^{wa,0). Then / = EZoMfH^'^'^ ^^ere d,(/) := 
(/, pI^'^^V Here and elsewhere in the following the convergence is in L^(wq;,/3)- 
Denote briefly 

(Lj * f)ix) := c„,^ y" ^ L^"'^(x, t)f{t)wc,0{t)dt, 

where ''^ are defined in (3.7). One can regard this as a nonstandard convolu- 
tion which is evidently associative but not commutative. We next show that the 
following Calderon type decomposition holds 

oo 

(3.13) f = Y^LfP*L';'^*f. 
Indeed, clearly for j > 1, 



Lr*f=i:^i^)Mf)Pi" 



:,/3) 

y=l 

and then 



L-/^*L-^*/ = 53a^(-;i^)d,(/)p(«./3). 

I/=l 

Consequently, 

oo oo 2^ 

^ Lf' * Lf' * / = d,{f)?t^<'^ + E E (^) '^^(/)pi"''^ 

j=0 j=l i/=l 

oo OO 

=rfo(/)pr) + EE«'(^)^^(/)p^"''^ 

oo 

=E^'^(/)p^^^ = /' 

where we used (3.6). Thus (3.13) is established. 

As we already mentioned L'^'^*L"'^*f = {L'^'^*L'^'^)*f. Since L"'^{x, t)L'^'^{t, y) 
is a polynomial of degree at most 2-'+^ — 1 in f and the quadrature formula (3.2) is 
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exact for such polynomials, we have 

Here we also used that L"''^{x,y) = L"''^{y,x) and the definition of the frame 
elements in (3.8). Consequently, 

which along with (3.13) yields (3.11). 

For the proof of (3.12) we denote Fjf := J2'j=o J2^eXj Evidently, 

(/,i^7/) = EEK/'V'£)r 

and passing to the limit as J ^ cxd we arrive at (3.12). □ 
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